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A Wilson loop is evaluated within the stochastic vacuum model for the case when the respective 
CNj ■ contour is self-intersecting and its size does not exceed the correlation length of the vacuum. The 

result has the form of a certain functional of the tensor area. It is similar to that for the non- 
t^j- ! self-intersecting loop only when the contour is a plane one. Even for such a contour, the obtained 

expression depends on the ratio of two functions parametrizing the bilocal field strength correlator 

taken at the origin, which is not so for the case of non-self-intersecting contour. 
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The small non-self-intersecting Wilson loop in QCD [], 

(W(C)) ee -L (trV exp (ig f dx^T^ 



has been evaluated long time ago. First, it followed indirectly from the respective quark- ant iquark 
potential emerging due to the dipole interaction with the external colourelectric field |1[ and then 
was found explicitly by performing the direct calculation in QCD ||. The term "small" here 
means that the typical size of the loop does not exceed the correlation length of QCD vacuum, T g , 
i.e. the distance at which the bilocal field strength correlator in stochastic vacuum model [0, |3| 
decreases. This length has been measured in the lattice experiments [[|, (see Refs. [7J for 
reviews) with the result T g ~ 0.34 fm for the realistic case of SU(3) full QCD. As a consequence of 
smallness of the loop w.r.t. T g , the field strength correlator of stochastic vacuum model can with 
a good accuracy be approximated during the evaluation of such a loop by gluonic condensate. 
Therefore for small loops the QCD vacuum is viewed as that of QCD sum rules |§ characterized 
by the infinite correlation length. However, self-intersecting loops are also of a great importance 
for QCD, since these are those loops at which loop equations [[| are nontrivial. In particular, in 
2D QCD such loops have been comprehensively studied in Ref. [^UJ (see Ref. [|11[] for a review). 
In the present letter, by combining stochastic vacuum model with the loop space approach we 
shall evaluate small self-intersecting Wilson loop in 4D QCD with arbitrary number of colours. 
As we shall eventually see, the resulting expression differs significantly from that for a non-self- 
intersecting loop. 

The idea we are going to employ is based on the possibility to represent the loop-space Lapla- 



cian 12 



1 ct+0 

A = / da I da' 



Sx^a'^x^a) 



standing on the L.H.S. of the loop equations in the following form 13 



A = i dx^a) v.p. [ da'x u {a')- f , (2) 

CT-0 1 

where v.p./ da 1 = J da' + J da'. It is worth mentioning once more that as it follows from 

CT+0 

the loop equations, the result of the action of the loop-space Laplacian onto the Wilson loop is 
nonvanishing only provided that this loop is self- intersecting. By virtue of Eq. (0), we get for such 
a loop the following equation: 

A (W(C)) = ~j dx^a) v.p. J da'x u {a') tr (F Xf ,(x)$(x, x')F uX (x')$(x' , x)) . (3) 

Here, A is defined by Eq. ([!]), x = x(a), x' = x(a'), F^ u = d^A u — d u A^ — iglA^, A v ] stands for 
the Yang-Mills field strength tensor, with A^ = A^T 11 , and $(x, y) = -^-V exp ^ig J A^iujdu^j is 

a parallel transporter factor along the respective part of the contour C . We can now take into 

1 From now on, we shall for brevity call by self-intersecting loop the loop defined at a self-intersecting contour 
and imply by the size of the loop the size of this contour. Besides that, all the investigations will be performed in 
the Euclidean space-time. 
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account the fact that the characteristic size of the countour C under consideration does not exceed 
T g 0. This enables us to use for the field strength correlator standing on the R.H.S. of Eq. (||) 
the expression known from the stochastic vacuum model || [|. In fact, for such a small contour 
joining the points x and x' the result for the bilocal field strength correlator suggested by this 
model is independent of the form of the contour and reads 

£(F fa ,(x)${x,x , )F Xp (x , )$( X r,x)) = ^f k {{6 ftX 5 vp -5 flp 8 vX )V{x-x')+ 



+- 



1 r 



(FMx - x')\5 vp - (x - x') p 5yx) + 8%{{x - x') p 5^x - (x - x') x 5 IX p) V x (x - x) 



(4) 



In what follows, we shall consider only nonperturbative parts of the functions D and Pi, which 
were measured on the lattice in Refs. [§, |], ||, [?]]. That is because the perturbative contributions to 
these functions, relevant to the UV divergencies ||14|| , yield a renormalization factor, which cancels 
out with the same factor appearing on the L.H.S. of Eq. (^]) during the direct renormalization of 
the Wilson loop Q. Therefore, from now on we shall deal with the renormalized Wilson loop as 
well as the renormalized charge g. Owing to Eq. (f|), Eq. (^) takes the form AVFfa;] = J[x], where 
W[x\ = (W(C)} and J[x\ = <f dx, 



F(x) 



^(cr) v.p. / da'x^(a')F(x — x') with 
1 



[6V(x) + 4Z>i(x) + x^Vxix)] . 



(5) 



The above equation for the Wilson loop can be solved by virtue of the method of inversion of 
the loop-space Laplacian proposed in Ref. [0 0. The idea of this method is to replace the original 

loop-space Laplacian (P by the smeared one = J da f da'G(a-a') 



s 2 



, , r, where G(a-a') 
is a certain smearing function. Such a smeared Laplacian can be inverted, which yields 



W[x\ 



dA (J 



x 



(G) 



J 



(G) 



(6) 



Here, 



<0[£]>f 



g(Q)=€(i) 

e(o)=s(i) 



is the Gaussian average over loops with the action 5* = \ J da J da'^{a)G 1 (a — a')£(a'), where 



G (a — a') denotes the inverse operator. Following Ref. fllql , we shall choose G(a — a') 
in the form G(a — a') = e~\ a ~ a '\l £ ) e \, after which the above action becomes local: 

S = \]da (ee{a) + \e{°)). 



More rigorously, this means that the area of the minimal surface spanned by C is not larger than Tg and that 
C does not contain appendix-shaped pieces of the vanishing area, but large length. 

3 The multiplicative renormalizability of self-intersecting Wilson loop has been proved in Ref. ]l5| |. 

4 This method has been successfully applied in Ref. Q to the solution of the Cauchy problem for the loop 
equation in turbulence. 
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Next, the averages on the R.H.S. of Eq. (||) are similar to those which one carries out during 
the reproduction of the one-gluon-exchange diagram within the same method when it is applied to 
the usual loop equation in the large- N c limit |16|]. In particular, this follows from the fact that the 



integral operator present in J[x] accidentally has the same form as the one standing on the R.H.S. 
of the large- N c loop equation. In the limit e — > 0, the only nonvanishing contribution appears in 
the first of the two averages on the R.H.S. of Eq. (BT) and reads 



1 1 



W[x] = l-^JdAjdaJ da' {I - G(a - a'^x^x^a'jx 





d 4 p 



(27T)' 



■ cxp 



Ap 2 (l-G(a-a'))+ip(x(a)-x(a')) F(p). (7) 



Here, the factor e -A P 2 (i-G(a-a')) is the result of the average (^ e i^p(^)-^'))^ ; anc i p(p) = 

J d A xe~ ipx F(x) is the Fourier image of the function F(x). Next, in Eq. (0), the factor (1— G(a— a')) 
in the preexponent was introduced in order to make out of the full integral over a' the principal- 
value one in the limit e — > 0. This factor disappears upon the introduction instead of A the new 
integration variable a = hr 2 + 2A(1 — G(cr — a')), where A stands for the UV momentum cutoff. 
Sending A to infinity we arrive at the following expression: 

W[x] = l--faf K f^— F(p) = 

l r , /•.,/•., F(z) 



dx^ j> dx'^ J d 4 z- 



8vr 2 J M J (z-(x-x')) 2 ' v ' 

The infinite integral (i.e. the integral over z or p) can be calculated in the small-loop case 
under study. In this limit, one can replace T>(x) and T>i(x) by their values at the origin, which 
according to Eq. (H) are related to each other as 



P(0)+2? 1 (0) = ^tr(i^(0)>. (9) 



On the other hand, according to the lattice measurements ]4], |5|, || [Tj, X'i(O) = aV(0), where 



a ~ 0.2 ±0.1 (see also Ref. [T3] for the discussion of this value of a). This yields the following 



approximate constant value of the function F(z): F(z) ~ l^fifjy^tr \-^L(0)) = The remain- 
ing infinite integral can easily be calculated in the limit T g 3> \x — x'\ under study e.g. from the 
first equality on the R.H.S. of Eq. (||). We have 



ipX 



-F(p) = C / d A P 5(p) 



d A p e 
(2tt) 4 ^2 

where A = x — x', and thus 

W\x] 



p 



16tt 2 



ds / d A pe 



ipX—p 



- ( ll. 



4A 2 



1-e T 9 



1 — i> dx,, & dx[, 



(x — x'Y 



1 - e 
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Expanding the exponential in this formula we finally arrive at the following leading nontrivial 
contribution to the Wilson loop: 

wm)^-t££t*M*>» (10) 

where = <f dx^x v is the tensor area corresponding to the contour C. Note that T g dropped 
out from this leading term, as it could be expected from the beginning, since l/T g was considered 
as an IR cutoff. 

The obtained expression for the Wilson loop is now worth to be compared with the respective 
expression for the small non- self-intersecting Wilson loop. In that case owing to the non-Abelian 
Stokes theorem and Eq. (f|) one has 



\ £min[' 



dxr^ix) J da Xp {x'){F llv {x)^{x : x')F Xp (x')^{x',x)) 

in [C] S min [C] 



/ 

exp 



9 2 



tr (-F^(O)) J da^ u (x) J da^(x' 



48iV c 

\ ^min[C] S m i n [C] 

Here, x = x(£) is the vector parametrizing the surface of the minimal area spanned by the contour 
C, S min [C], and £ = (£\£ 2 ) stands for the 2D-coordinate. Next, do^ u (x) = \fg(C)t flu (t;)d 2 t;, 
where #(£) is the determinant of the induced metric tensor g ab (C,) = {d a x ll (^))(d b x l j i (C,)) and 

tfj,v(0 = £ab (da x /j.(0) {9b x v(0) I \] 9{£) * s ^ ne extrinsic curvature tensor. Since the contour C under 
discussion is very small (and consequently the same is E min [C]), the points x and x' are located 
very closely to each other, and therefore ^(0^w(O — tfiuiQ = 2- Finally, taking into account 
that/d 2 £ v /^(0= Areaof£ min [C] = S min , we obtain @ (W{C)) ~l--^tr (F%(0)) S^. For 
a plane contour, this expression has the form similar to our Eq. (|10D, since for such a contour 

In this case, the main difference between these two expressions stems from the a-dependence of 
Eq. (|T0|). This dependence is due to the fact that the functions V and V 1 contribute to the self- 
intersecting loop in the nontrivial combination (|5D, whereas the non-self-intersecting loop depends 
only on their sum at the origin, expressible owing to Eq. (^) via the condensate alone. Another 
obvious difference of the two expressions for Wilson loops is that the tensor area for self- intersecting 
countour can be vanishingly small even for a very large contour (although we do not consider such 
contours) and even vanish completely for the eight-shaped contour with equal petals, whereas for 
a non-self-intersecting contour it could vanish only together with the contour itself. Moreover, 
in this respect it is worth pointing out once more that the comparison of the results for self- 
intersecting and non-self-intersecting Wilson loops is only possible for plane contours, since for 
non-plane ones Eq. (pTTj) is not valid. 

In conclusion, by making use of the method of inversion of the loop-space Laplacian, we have 
restored a small self-intersecting Wilson loop from the bilocal field strength correlator of stochastic 
vacuum model. There turned out to be two main differences of the obtained result ([!(]) from that 
of the non-self-intersecting loop. Firstly, Eq. ( [10 ) depends on the tensor area of the contour, 
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rather than on the area of the minimal surface, and can therefore vanish for some class of contours 
{e.g. for plane eight-shaped contours with equal petals). Secondly, the obtained result depends 
on the ratio of nonperturbative parts of the functions T>i(x) and V(x) at the origin, which is not 
the case for a non-self-intersecting contour. However, for plane contours, the functional form of 
the obtained result coincides with that of a small non-self-intersecting loop when the latter one is 
expressed in terms of the tensor area. 
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